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Wheeler-DeWitt approach and ordinary QFT. Remarkably, the Penrose-Hawking singularity of usual Friedman- 
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and compared with loop quantum cosmology. It is argued that it is the regularization of the classical Hamiltonian 
performed in this last theory what avoid the singularity, rather than quantum effects as in our case. 
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In accordance with the celebrated Penrose-Hawking singularity theorem, Friedman-Robertson-Walker (FRW) cosmologies 
give rise to a singularity when the strong energy condition holds, p + 3p > 0, p being energy density and ^pressure flJJ] . A simple 
way to avoid this singularity is introducing a scalar field that breaks the strong energy condition [2, 3, 4]. Another possibility 
0^ 1 is to considerquantum effects due to vacuum polarization, as the one due to a massless scalar field conformally coupled with 
gravity (5|, la |7[] . Remarkably, there is a very natural and fundamental alternative: to depart from the Wheeler-DeWitt equation, 
HQ — HIH], with H the quantum Hamiltonian, by assuming that time has an absolute meaning. The question of a possible 
singularity is then addressed in terms of a Schrodinger equation with additional conditions, what specifically defines our theory: 

CT 

u' ihdMt) = H$(t), $(**) = *, (H)* = o, ||*|| = i. (i) 

bJO, 

We thus plainly assume the choice of a time direction to have a physical solution, what goes against common lore. In this sense 
I [ our approach is absolutely original, revolutionary and, as we will prove, clarifying and predictive. 
■ Some basic technical details on what will be demonstrated below. The quantum Hamiltonian H, obtained with the usual rules 
of quantum mechanics, is generically symmetric but not self-adjoint. By von Neumann's theorem ifiolflltl . it can be extended to 
a self-adjoint operator (sometimes in infinitely many ways). Stone's theorem then applies, leading to a solution valid at all time t 
and, consequently, we can compute the average of the quantum operator a corresponding to the classical scale factor a. That is, 
we compute the following effective scale factor a e //(t) = ($(t)|a$(t)), where $(t) is the solution of the effective Schrodinger 
t-H ■ equation above. It is not difficult to see that if $(<) belongs to the domain of the operator a at any time, then the effective scale 
factor is always strictly positive, and we can conclude that the singularity is avoided. Physically, the self-adjoint extension of the 
Hamiltonian operators that appears in FRW cosmologies can be understood assuming that there is an infinitely barrier potential 
at the point a = 0, then when the effective factor scales approaches to zero, at some finite time, it bounces and grows. 

However to compute averages one usually works in Heisenberg's picture, so A = i [H, A], for A any operator involved in the 
calculation. But, using this formula one turns out to obtain, at some finite time, a negative value for the average of the scalar 
factor operator. Such a contradictory result can be explained by the fact that, at some finite time, the commutator between the 
Hamiltonian and the operator A is not well defined, what invalidates the final result (physically one can explain this taking into 
account that in Heisenberg picture the boundary conditions do not appear, i.e., the barrier potential are not introduced and then 
the effective factor scale has the freedom to take all the values in R). That is, the average of the scale factor is positive, but we 
do not have any method to obtain an analytic information about its behavior, because the Heisenberg picture fails to work, and 
it also turns out to be impossible to obtain an explicit solution of the effective Schrodinger equation. To this end loop quantum 
cosmology (LQC) will be invoked lfl2l[T3Tl . In what follows, we present a simple demonstration of the above approach and will 
explicitly see how this theory avoids the singularity. It will be also shown that it is the regularization of the classical Hamiltonian 
iffl [T5I llcjll what avoid the singularity, rather than quantum effects. 

In the first of the three Appendix in the paper we present a brief mathematical review about the theory of self-adjoint extensions 
of symmetric operators based on Von Neumann's theorem. In the second one, we apply the effective formulation to the case 
of a barotropic fluid where on can see clearly the physical meaning of the self-adjoint extension of a symmetric operator. As 
specific examples, the dust and radiation cases are treated in detail, showing that the self-adjoint extensions of the respective 
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Hamiltonian operators can be understood assuming that there is an infinitely barrier of potential at a = 0. Finally in the last one, 
we show (resp. review) how to derive the standard quantum fields theory in curved space-time from the effective formulation 
(resp. the Wheeler-DeWitt equation). We also obtain, from the effective formulation, the semi-classical Einstein equation, that 
is, the back-reaction equation. 

n. THE PROBLEM 

In this Section we consider an homogeneous and isotropic gravitational field minimally coupled to an homogeneous scalar 
field, which Lagrangian is given by J3l 

L{t) = l^ (c2fc " k2)a + ^ V " F(0)a3 ' (2) 

where G is Newton's constant and k the three-dimensional curvature. We are interested in the case k = and V = 0, previously 
studied in lfl7i [Till within the framework of LQC. This interest comes from the fact that in the chaotic inflationary model with 
V = \m 2 cf) 2 , at very early times before the inflationary period, one has <fi ^> V (see for details 11910 . Then the potential can be 
neglected and one has p = p. Consequently this model give rise to a singularity at very early times, that we want to avoid using 
the effective formulation described in the Introduction. 

Defining the angle variable ip by <f> = y/3c 2 / AttG ip, the Lagrangian becomes (l p denotes the Planck length) 

Using the conjugate momenta, p a = —j 2 aa andp^ = j 2 if;a 3 , we can write the Hamiltonian as 



H =^[-(aPa) 2 + P l]. (4) 



The classical dynamic equations are 



1 ■ . -1 

ap a , (apa) = 3H, ip = —zr-zPii, Pi> = 0, (5) 



together with the constraint H = 0, that is (ap a ) 2 = p 2 ,- Integrating ((5]) we obtain the following solution 

> a(t)p a (t) = p* a , p^±{t) EE pj ± = ±p* a , (6) 
a(t) = a* 



1 _3p* a (t-t*) 



7 2 (a*) s 



1/3 /.N 

a 



Note that the solution (a(t),ip±(t)) is defined in the interval (t s , +00), where £ s = t* + {"f 2 /3p* a ) (a*) 3 . At this time 
we have a(t s ) = 0, and 4>±{t 8 ) — ±00, that is, the dynamics is singular at t — t s . Note that we can write a(t) = 
a* [1 — (t — t*)/(t s — t*)] 1 ^ 3 . Finally, we see that from Eqs. Q we have dlaa/d^± — =Fl, and conclude: 

o = a*e =F W ± -^. (7) 
A. Quantum dynamics 



We now use the quantization rule: 



9 AB PAPB — -h 2 V A V A = —^=d A (^\g~\g AB d B ), (8) 



1.91 



and obtain the quantum Hamiltonian 



H=^-^(ad a ad a -8fc). (9) 
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Introducing the operators ap~a = —iha 1 / 2 d a a 3 / 2 ] p^ = —ihd^,, we can write H = -^a 3 ^ 2 [p 2 b — (ap a ) 2 ]a 3//2 . The 
dynamical equations in the Heisenberg picture are 

dip 1 » dp^ 

with ($|*) = J °° da J R d^*(a,ip)a 2 ^(a,ijj) as inner product and operator average (A)^ = (^\A^), ||^|| = 1. 

Example 1. Consider the wave-function (|p*| = \p*^\): V(a,ip) = e~ ^^'^ e~ <dL ^r L - e i('"( a / a K+^;). Then: 

(a)w = ae CT / 4 = a*, = ip*, (ap^)* = p*, (p^>* = and = 0. 

B. The Wheeler-DeWitt equation 

Compare at this point with the Wheeler-DeWitt equation paradigm (WDW) JH] 

ff$ = 0^ (ad a ad a -d 2 2 )<S> = 0, (11) 
with general solution (a is a length constant) 

*(o, </>) = f+(Ha/a) + + /_(ln(a/3) - V)- (12) 

The quantum version of Eq. |7]) is $+(0, -0) = f + (ln(a/a*)+ip — ip*), (a, ^) = /_ (ln(a/a*)~ijj + ip*), with /± a function 
picked around 0, as for instance f±(z) = e~ z 1° ' . From this result we see that the wave is always picked around the classical 
solution, but we cannot conclude that its dynamical behavior is singular since here time does not appear. 

In order to understand the dynamics, we postulate the effective equation (HJ. Note that (H)^ = implies (H)&u-\ — 0, 
and 1 1^| I = 1 implies ||$(t)|| = 1, Vf G M. Then, if the solution of the problem exists for all t, it is easy to prove that 
the effective scalar factor, a e ff(t) = (&W t ), never vanishes. In fact, the condition \\^\\ = 1 implies ||$(f)|| = 1, i.e., 
J °° da f R dijja 2 \<S>(t, a,ip)\ 2 = 1, and thus we have: a e //(t) = (a)$(t) = / °° da J R dipa 3 \<&(t, a, tp)\ 2 =/= 0. 

To do the calculation, we consider the quantity (a 3 )$( t ). We have ^(a 3 )*(t) = a 3 ])$(t) = — ^(oPa)*(t) and 

^t 7 (d 3 )<s>(t) = — y^{[H, 5p^])$(t) = — T-(-ff)$(t) = 0, due to the remark above. Consequently, we obtain 

(a 3 )$ (t) = (a 3 )* - ^<^>*(*-**)- (13) 



For the function of Example 1 , we get 



<a 3 ) $(t) =a 3 e 9 -/ 4 -^K(t-r), (14) 



and this contradicts the fact that (a 3 )$( t ) > 0,Vt G K. However, since the operator H is symmetric and real, using von 
Neumann's theorem 030] it can be extended to a self-adjoint one, and then the solution of the problem (Q~|) exists here for any 
t (Stone's theorem). From this result, we conclude that there is a value of t for which some of the commutators [H , a 3 ] and/or 
[H, apa] do not exist; thus the final result ( TPfl ) is incorrect. The drawback of this method is the lack of an analytic procedure to 
calculate the average since, in general, there is no explicit formula that gives information on the regular behavior of the average 
of the scale factor operator. Fortunately, a useful way exists to directly analyze the singularity, namely loop quantum cosmology 
(LQC). Before using it in our problem, we consider another example where the above contradictions can be easily depicted. 

Example 2. Consider now the problem 

iflt$ = -ihcd x ® = cp$, y x G [0, 2tt], $(0) = V, (15) 



being p self-adjoint in the domain [11]: Dp — {ty absolutely continuous in [0, 2ir], d x ^$> G £ 2 [0,27r], 

^f(0) = 4"(27r)}. Let $(t) be the solution of our effective formulation (fTBI l. We want to calculate (x)^m = J„ 27r x\$>(t, x)\ 2 dx. 
Using \p, x] = —ih, we get {x)$i t \ = c, i.e., (x)$m = (x)s, + ct which is not positive Vi. What actually happens is that, for 



4 



some t, we have x$(t) ^ Dp, and then px$>(t) has no sense, neither the formula (i)$(t) = i([cp, x]Wt). To see this in detail, 
consider the initial state 



, - i x, for a; € fO, 7rl, ,, „ 

* (x) = ^ for , 6 i 4 (16) 



Fourier analysis provides the following solution of the Schrodinger equation 



r._4 F ^_ 

I 2 7r ,ife( 2n+1 ) 



Then, at £ = we have x$(0, x) S Dp but if we choose t = n/c we obtain x<I>(7r/c, 5c)l_ = 0, and x$(7r/c, x)| x=27r = \/6tt/c, 
what means effectively that x§(ir/c, x) ^ Dp. However, note that (x)&m exists for all t G M, its value being: 

o < (*>•(*) = y o ^{f--E(^Ti? C0S[(2n+1)(a; - rf)] ) dx<2 "- 



III. LOOP QUANTUM COSMOLOGY TO RESCUE 

We shall now involve (a simplified version of) LQC (for a rigorous formulation see lfl3l [l6lo . with different variables and a 
different quantum space of states, adapted to make contact with our theory above. 

Consider the variables p = a 2 and x = a. Their Poisson bracket is {x,p} — = \. We also consider the holonomies 

hj(n) = e~ % ^ aj — cos(^) — ictj sin( 2 ^_) l2lll . where <Xj are the Pauli matrices and i is the Barbero-Immirzi parameter. We 
easily obtain Ashtekar-Barbero's formula B21I1 



a ^ = TwTrT.vMWj 1 ®,*}- (17) 
W V j=i 

To get the gravitational part of the Hamiltonian, we cannot directly use this one: H grav — — §^qX 2 ^/p, which leads to singular 
classical dynamics. We may use the general formulae of loop quantum gravity (LQG) to obtain the regularized Hamiltonian 

*W = ~3^W E £ijkTr Ml)h j (l)hr\l)hj\l)h k (l){h^(l),a 3 }] = -^sin 2 ^, 

p i,j,k 

which is bounded when the extrinsic curvature x/2 (a half of the velocity of the scalar factor) diverges, and approaches H grav 
for small values of x. Then, taking this regularized Hamiltonian as the gravitational part of the full one, this last is given by 

Si 

'V C o LX 1 q o 

Hl = ~ 2^ ashl ~ + 2^ a p *> (18) 

and the dynamical equations are 

c 2lx 2 

a = {a, H L } = — sin , x = {x, H L } = t^pI- (19) 

Imposing the Hamiltonian constraint H t = 0, we obtain 



and since p^ = p*^ is constant, we get the following bounce, d = when a = -y = 2l p y Consequently, there is no 
singularity because the range of a(t) is 



2l p \j -jf-, +oo ) , Vt € M. In fact, at earlier times the scalar factor is very big, then it 
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decreases, and when it arrives at the turning value it increases forever. Moreover, this solution yields a period of inflation [22], 
namely, from the Friedmann Eq. (l20l i: a > 0, for a £ ^2l p \J~^jf- , 21 p \j —^^j ■ Finally, note that when a ^> l p , Eq. ( f20b 
coincides with ©. 

A different way to understand these features is to write Eq. (|20|) as (-) = -%p e //, where we have introduced the effective 

2 / 2 2 \ 

energy density p e ff = ( 1 ifr^ ) • Taking the derivative, p e ff — —3 (~) (p e ff + Peff), where the effective pressure 



is p e f f = -^p ( 1 — 3 7 4^4 e 2 j • But then is easy to see that the strong energy condition p e f f + 3p e f / > is broken, when 
the scale factor lies in the interval f 2l p \J^^- , 2l p \J —^^-J , consequently, the singularity is avoided. Moreover, for a € 



^ pV /^p, 2l p y ~$~^f~J ' th ere is a period of super-inflation; that is, in this interval, one has yjj < — 1. 

The following remark is in order. Similar results are obtained in the case k — 1 and V = 0. Now the classical Hamiltonian is 
given by 

H=^- s (-(ap a ) 2 +pl- 1 4 c 2 a 4 ), (21) 



and the regularized one, is 



2" 

Then, using the Hamiltonian constraint H L = 0, is easy to obtain the Friedmann equation 



" 2 =[^- cZ )[ 1 + ^^)> (23) 



2 \ / 2 2 

Pip 2 \ / i , 2 Pi/) 1 



and since, p^, = p\ is constant and a 2 > 0, we can deduce that a E 



and clearly the singularity is 



avoided. In this case, we have an oscillating universe. 

Another equivalent way to do this, is to use the variable x = a + c, then following rf25[ l26ll we obtain the regularized 
hamiltonian 

H t ee -t^a (sin 2 (^^) ~ sin 2 ( t ) + 2, 2 ) + -^a" 3 ^ . (24) 
It is clear, that from this last regularized hamiltonian the scale factor has the same behavior that described in equation d23V 



A. Quantization 



To quantize we perform the usual change {A, B} — » —jr[A,B}. Note that the system is ^^-periodic with re- 
spect to the variable x, thus we consider the space of : ^-periodic functions and introduce the inner-product ( , 3 / | < E > ) ee 

2-kc 

J_ dip f 2„ c dx^*(x, iji)$(x, ip). Completion of this space with respect to this product is the space of square-integrable func- 
tions in [— ^t 2 , 2 -f £ ]- Note that, rigorously, the definition of the Hilbert space is more complicated: C 2 {M.Bohn d^Bohr) where 
^Bohr is the compactification of R and fiBohr the Haar measure on it 11811 . However, for our purposes the Hilbert space 
£ 2 [-^, will suffice. 

1/2 



We quantize the variable p as above and, using the fact that p > 0, we can define p 
V ee p 3 / 2 , and the scale factor a = p 1 / 2 . The eigenfunctions of these operators are \n) ; 

(p) n = fc\n\l 2 p , (V) n = (f i\n\l 2 p f' 2 and (a) n = Jfi\n\l p . Using Eq. O, 



, the volume operator 



— e 2<= x , and their eigenvalues are 



47TC 



^l 2 l 



Trj2oMl)[h-Hl),V l, % 

3=1 



(25) 
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and consequently the corresponding quantum operator is 




(26) 



whose eigenvalues, when n 3> 1, satisfy (a ) n = l/(a) n . 

The quantization of the gravitational part of the Hamiltonian, depends on the order we fix. For instance, H gr av.i 



:52ji 



^wE^fc^r h-\l)h-\l)h k {l)[h- k \l),V]Hl) fej(l) 



gives us a self-adjoint operator, or the direct quantiza- 



tion of the expression — ^§-asin 2 — yields 
If we use this operator d27l) as the gravitational part of the full Hamiltonian, then this is given by 



//,„,„.,: --^a^sin 2 ^)^/ 2 . (27) 



- -^f« 1/2sin2 (7) « V2 + ^ ^fv% (28) 

and in this case the WDW equation becomes H t <I> = which, expanding $ as $ — EneN ^(VOI 7 " 1 )' turns into 

2^$„ - |n(n - ^l 1 /**^ - | n (n + ^l 1 / 4 ^ 

+4 + l\-y/\ n - iff 9 2 2 $„ = 0, neN. (29) 

Summing up, the effective equation ihdt& = H L <& with the condition (H L )$,( t ) = yields an average of the scalar factor 
operator that has essentially the same behavior as the classical solution of Eq. ( 1201 ). This owes to the fact that the domain of the 
holonomy operators is the whole space, so that one can safely use the Heisenberg picture in order to obtain the quantum version 
of the classical equations. This gives generically small corrections to the classical behavior. 

A final remark is in order. The singularity is avoided in the classical theory after regularization of the Hamiltonian. Quanti- 
zation of this new Hamiltonian provides then a self-adjoint operator. It is important to realize that it is the regularization of the 
classical Hamiltonian what avoids the singularity, rather than the quantum effects. This is overlooked in some papers, where it is 
claimed that quantum effects are essential to avoid the big bang singularity Note that in these approximations one 

already starts from the quantum theory and then, using the quantum operators an effective Hamiltonian is obtained ll27i EH E9I1 
which, in fact, is in essence the Hamiltonian ( TT8l . This is maybe the reason why it is plainly concluded there that quantum 
effects, provided by LQC, are responsible for avoiding the big bang singularity. Here, with our alternative formulation we have 
shown, by means of explicit examples, that this need not be the case. 



IV. CONCLUSIONS 



We have presented here an effective formulation that naturally avoids the big bang singularity: in essence Schrodinger's 
equation with the condition that the average of the Hamiltonian operator be zero. This is different from the Wheeler-DeWitt 
equation where one impose that the Hamiltonian operator annihilates the wave-function, and the arrow of time is yet to be 
selected. In our theory, physical time has essentially the same meaning as in the classical theory, and the relevant quantities 
are averages of quantum operators, as e.g. the average of the scale factor operator — which is by definition strictly positive — 
and no singularity appears at finite time. Our approach is remarkably natural (once time is assumed to exist), revolutionary 
and predictive, albeit rather non-trivial. It does not seem easy to produce an analytic formula that provides information on the 
behavior of the observable averages. Only numerical results look feasible at this point. 

Another way to deal with the classical big bang singularity is LQC. We have here involved a simplified version of this theory 
and shown that, in contradistinction with the theory presented above, in LQC it is the regularization of the classical Hamiltonian 
that seems to avoid the singularity, and not the quantum effects obtained after quantization of the regularized Hamiltonian. 



V. APPENDIX A: SELF- AD JOINT EXTENSIONS OF SYMMETRIC OPERATORS 



In this mathematical Appendix we present a brief review of the theory of the self-adjoint extensions of symmetric operators. 
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Let A be a linear operator that is defined on a dense subset of a separable Hilbert space TL . The adjoint A* of A is defined 
on those vectors $ S for which there exist $ G such that ($| A*) = ($|*) V* e Z)^, and A* is defined on such $ as 

The graph of an operator A is a subset of 7i (J) 7Y, defined by = {(<!>, A$); $ S -D^}, and A is called closed, which 

is written as A = A, if its graph is a closed set. An extension of an operator A, namely A ext , is an operator that satisfies 

D A c D A ext and = A ® V$ e D A- 

An operator A is symmetric if ($| A^) = (A$|'I') V^P, $ £ -D^. Then, a symmetric operator A always admits a closure (a 

minimal closed extension), which is its double adjoint, i.e., A = A". The adjoint of a symmetric operator A is always a closed 
extension of it, and it is self-adjoint when = D i t . The deficiency subspaces Af± of the operator A are defined by 

A4 = {$€D At , A t $ = z±$, ±Im(z±)>o}, (30) 

and the deficiency indices n± of A are its dimensions. Note that, these two definitions do not depend on the values of z±. 
The following theorem is due to Von Neumann: 

For a closed symmetric operator A with deficiency indices n± there are three possibilities: 

a) If n + = Ti- = 0, then A is self-adjoint. 

b) If n + = ri- = n > 1, then A has infinitely many self-adjoint extensions parametrized by an unitary n x n matrix. Each 
unitary matrix U n : A/+ — > A/1, characterizes a self-adjoint extension Ajj n as the restriction of A^ to the domain 

^i Un = + +tfn$* + ;$e r> A $ z+ eA/+}. 

c) If n + 7^ then A has no self-adjoint extensions. 



VI. APPENDIX B: EFFECTIVE FORMULATION FOR A BAROTROPIC PERFECT FLUID 



In this Appendix we apply our effective formulation to the case of a barotropic perfect fluid with state equation p = up. The 
Lagrangian of the system in the flat case (k = 0) is 

i 2 

L = — ^-d 2 a - p(a)a 3 . (31) 

The momentum and the Hamiltonian are respectively p a = —j 2 da, and H — —-^r^Pa + p(a)a 3 - Using the conservation 
equation a = —3H(p + p) we have p(a) — po (a/ao) 3 ^ +1 \ then the dynamical equations become 

2 

a = 1-; p a = - § 2 + 3up(a)a 2 , (32) 

with the constraint H = 0. 

The quantization rule dHJ give us the following Hamiltonian operator 

H=^- r d 2 + p(a)a 3 , (33) 

which is symmetric with respect to the inner product (<I > | , I') = / °° daaQ* (a) 1 ^ (a) of the Hilbert space £ 2 ((0, oo), ada). 

To apply the theory presented in the Appendix A, first we consider the case u> = (dust matter), whose Hamiltonian is 
H = ^-d 2 a + po«o- To study the self-adjoint extensions of this operator we need to determine the deficiency subspaces A/±, 
that is, we must solve the equation HQ = z±$ with 1 1$| | < oo. Since the definition of these spaces do not depend on z±, we 
choose, z± = ±ipoa,Q. Then the solutions of H$ = ±ipoa,Q& are the Airy's functions $i,± = Ai(f3±a) and $2.± = Bi((3±a), 

2 3 1 /3 

where j3± = (^-j^r^J e ±!7r / 4 ll30Tl . However, only $i,± has finite norm, and then both spaces has dimension 1. Von 

Neumann's theorem says us that H has infinitely many self-adjoint extensions, namely Hsa, parametrized by an unitary lxl 
matrix, i.e., by e a being a 6 ML To obtain an explicit expression of the domain of these self-adjoint extensions we must impose 
USUI (H SA (<£> + + e M $_)|*) = ($+ + e la <5>-\H SA q)) V* E D &sa . It is not difficult to show that, this condition is 
accomplished when 

*(0) Ai(0) 1 

^ ^(0p + |l + tan(a/2r r ' Whh ? * ^ (34) 



8 



That is, for different values of r we obtain different self-adjoint extensions. Here a very natural extension is obtained choosing 
r = 0, that is, imposing \&(0) = 0. Physically, this is equivalent to assume that at a = there is a infinite potential barrier (in the 
same way that for no-relativistic one-dimensional barrier problems), then the existence of a solution all the time is guarantied 
because when the scale factor decreases to zero, at some finite time, the potential barrier forces it to grow. Moreover, this 
assumption explains why the Heisenberg picture fails to work, because in the Heisenberg picture the boundary conditions do 
not appear, and the effective scale factor has the freedom to take all the values in R, in particular, or negative values. We can 
conclude that if we want to work in Heisenberg picture we must introduce some kind of potential barriers that prevent that the 
effective scalar factor takes negative values. 

Once we have obtained a self-adjoint extension we apply the effective formulation ([TJ to the problem 

ihdMt) = JL-%$(t) + p a 3 m, (35) 

with the additional conditions $(£*) = 'J, (Hsa)^ =0, \\^\\ = 1, that gives us an strongly continuous unitary one- 
parameter group defined on C 2 ((0, oo), ado) (Stone's theorem), namely e~^ HsAt . The solution of our problem can be written 
as $(t) = e~^ HsA ^~ e ^ for all * e -D# SA satisfying (H S a)^ = and \\^\\ = 1. As an example of initial condition, if 
r = 0, one can take 

\E(a) = , a - e - h ^ 1 e^ lnia/!i)p '\ with p* = -J 2p {a af^e-l" - 7i 2 (25/4+ l/(2a)). (36) 
(an) / v 

For this initial state, the effective scale factor a e //(i) = (a)$(t) grows forever for t > t* in the similar way to the classical one 
(the classical limit holds far of the turning point a = 0). For t < t* the effective scale factor decreases to zero, but at some finite 
time it bounces, due to the potential barrier, and then it grows to infinity. 

Finally, we study the case u = 1/3 (radiation). The Hamiltonian is H = j^r^d 2 + an d me solutions of the 

equation HQ = ztipoa^Q are the Airy's functions Qi,± = Ai(f3±(a =p icta)) and <1>2,± = Bi((3±(a =p icto)), where 

(n 2 3 \ 1/3 

In this case the dimension of both deficiency subspaces is 1, then as the dust matter case, H has infinitely many self-adjoint 
extensions parametrized by an unitary lxl matrix, and the self-adjoint extensions are determined, once again, by the boundary 
condition ^(0) = r^'(0), with r?l. Now an initial condition for our effective formulation, that exhibits the same behavior as 
above for the effective scale factor, is given by the function 

^(a)= - a - e ~ 1 4^ £1 e i0"(°/a)p , ) i with p* = - J2p (af.a) 2 j 2 e-^ - ft 2 (25/4 + l/(2a)). (37) 

((TTt) 1 ' 4 v 

We finish this Appendix with the following remark. When the three-dimensional curvature is positive (k — 1), the Hamil- 
tonian of the system is H = —^aPa + p( a ) a3 ~ \'Y 2 c 2 a. Then the Hamiltonian constraint restricts the value of the scalar 

factor into the interval (0, A) with A = ^2p ao^ +1 V( c 7) 2 ) 3 " +1 > an d tn i s sa Y us tnat we must ta ke a s Hilbert space, the space 

C? ((0, A), ado). Now for oj < 1, a = is a regular singular point of the ordinary differential equation HQ = z±Q, then ap- 
plying the Frobenius method we can deduce that there exist two independent solutions of the differential equation, consequently 
both deficiency indices are 2, because the domain (0, A) is finite (excepts for lu = —1/3). Then the self-adjoint extensions are 
parametrized by an unitary 2x2 matrix, and the more natural boundary condition is to assume that the wave-functions vanish 
at two boundary points. Physically this means that the scale factor is confined in a very deep well potential, and we have an 
oscillating universe whose effective scalar factor never vanishes. 



VII. APPENDIX C: QFT IN CURVED SPACE-TIME FROM THE EFFECTIVE FORMULATION 



For the flat FRW universe, the action that describes a massive scalar field conformally coupled with gravity in the presence of 
a barotropic fluid , is given by 



S 



dt 



dx 



7 2 -2 / / \-3(w+l) 3 

—a a- pn (a/a a ) a 



a 6 £ d 



(38) 



with C 4> = ^ 



IKca 1 



(V</>) 2 



m_c i2 



uhc 3 R 2 ^ 2 where R — -^(d 2 + a'd) is the scalar curvature. (Note that in this 



Appendix <fi has energy units). Integrating with respect x and expanding <fr in Fourier series (cj) = J^k 



feez 3 



2TTi± 



") one obtains, 
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S = j M L(t)dt, with 



L(t) = L 3 { --La 2 a - p (a/ao)" 3(w+1) a 3 + a 3 £ }> , 

feez 3 



(39) 



where ^ = ^| - " - db^I 



Using now the conformal time dr\ = ^f-dt, (t p being the Planck time) and defining the function tp^ = \J -jjf-^4>j: we obtain 

47T 2 |fc| 2 



L(t)dt = ^-L{ri)dr], with 



L(„) 



(7) -*<«/«.>-+" 



(40) 



where we have introduced the Compton wavelenght l c = and we have defined po = ifPo- ^ ne important remark should be 
made that in this Lagrangian we have suppressed the terms — ^ ^V'g 

ft a' 



The conjugate momenta are p a = —f-^r, Pip % — ip'r, and the Hamiltonian is given by 



(41) 



where m p is the Planck mass and 



4 1 

t t / \ ~ii \-3(w+l) a 2/ \ 1 

17(a) = po (a/ ao) — , w g (a) = -5 



47T 2 |fc | 2 

L 2 



+ 



The quantum theory is obtained making the replacement p a — > — ihd a and p^^ — ► —ihd^ n . Then the quantum hamiltonian 
is given by 



5- " 2 - 



2m, 



^ 2 +f7(a) + ff m (a,V), 



(42) 



where the matter hamiltonian is i7 m (a, tp) = J^kez 3 (^S^g^fe + l^i) > wnere we nave introduced the creation and anihila- 
tion operators 



it EE —L_ 



(43) 



Now, we show how one can obtain the QFT in curved space-time from the WDW equation. If we consider the matter field as 
a small perturbation, we look for solutions of the WDW equation with the form $(a, tp) = ^(a)x(a, tp). After substitution in 
the WDW equation we obtain: 



2m. 



-d 2 2 * + [/(a)# 



X + 



*^%>X + — d a ^d a X + ^H mX 

2m p m p 



= 0. 



We assume at this point that * is the solution of the equation 

ft 2 



2m, 



-U{a)^ = 0, 



(44) 



(45) 



and we make the change ^ = e r s , then we obtain the system 



^-m + ftfts = 



(46) 
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To solve this equations we neglect, as Rubakov does 0211 . the second derivative with respect to a, then we obtain the system 

[d a S) 

(47) 



M_c/ (a) = o 

-ih^OaX + HrnX = 0. 



The first equation is the classical Hamilton-Jacobi equation, and the second one is the quantum Schrodinger equation that can 
be solved choosing as solution of the Hamilton-Jacobi equation S(a) = y / 2m p U(a)da, and introducing the conformal time 
^ = then the Schrodinger equation becomes ihd T x = H m {a{r), ip)x- 

Finally, we device a method to obtain the QFT in curved space-time from the effective equation ihd v § = HQ. Assuming that 
the matter field is an small perturbation, we look for solutions of the form <&(a, ip; rj) — ^(a; r])x(ip', rj) where ^ is the solution 
of the equation 

h 2 

ihd n V = 7^~d 2 a^ + U{a)y, (48) 

and we assume that ^ is a function concentrated around a classical solution, namely a c (rj), of the following equation 

1 



2m p 



pi + U(a) = 0. (49) 



By inserting $ in the effective equation one obtains ^>ihd v x = ^H m {a, ip)x> an d since ^ is concentrated around the classical 
solution, one can approximate ^>H m (a, tp) by ^H m (a c (r]), ip), and then one obtains ihd v x — H m (a c (r]),ip)x- 

We end with a last remark. From the effective formulation it's not difficult to obtain the semi-classical Einstein equations. 

Effectively, starting with the condition (H)<s> = 0, if we take the wave function used above (now picked around a c + Sa c ), one 
approximately obtain 

Pl c +Sa c + U{a c + 8a c ) + (H m {a c (r)) + 5a c {rj), ip)) x ,ren = 0, (50) 



2 m p 



where the quantity (H m (a c {rj) + 8a c (rj),tp)) x has been renormalized. 

Since a c is solution of the equation ( |49l , one also obtains, in the linear approximation, the following back-reaction equation: 



-^-a' c (Sa c )' + U'(a c )Sa c + {H m (a c (rj), tp)) x , ren = 0. (51) 



Finally, observe that the derivation of the semi-classical Einstein equation from the WDW one is not a completely clear case (see 
for example lO ). 
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